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NONLINEAR  OSCILLATIONS  IN  EQUATIONS  WITH  DELAYS 

Jack  K.  Hale 

1.  Introduction.  These  lectures  are  concerned  only  with  some 
aspects  of  bifurcation  theory  in  the  local  theory  of  nonlinear 
oscillations  in  equations  with  delays ; that  is,  behavior  of  solu- 
tions near  an  equilibrium.  In  particular,  we  study  how  the 
qualitative  behavior  of  solutions  change  as  parameters  vary. 

A detailed  study  of  the  local  theory  is  important  in  order  to 
know  the  types  of  solutions  to  expect  in  a global  problem.  Of 
course,  there  is  no  reason  to  only  study  local  theory  near  an 
equilibrium.  One  should  study  how  the  qualitative  behavior 
changes  near  any  invariant  set  - for  example,  behavior  near  a 
periodic  orbit,  behavior  near  an  orbit  which  connects  a 
saddle  point  to  itself,  etc.  More  complicated  behavior  is  ex- 
pected near  these  large  invariant  sets.  One  can  obtain  in- 
variant torii,  homoclines  points  which  exhibit  a chaotic 
behavior,  etc.  Ke  restrict  ourselves  in  these  lectures  to 
behavior  near  equilibrium. 

The  simplest  type  of  smooth  bifurcation  is  from  an 
equilibrium  to  a periodic  orbit  - the  so-called  llopf  bifurca- 
tion. In  Section  2,  we  discuss  the  Hopf  bifurcation  in 
equations  with  finite  delays  permitting  the  bifurcation 
parameters  to  be  the  delays  themselves.  At  first  glance, 


) 


such  a result  does  not  seem  possible  because  the  vector  field 
in  the  equation  is  not  differentiable  in  the  parameters.  The 


theorem  does  require  some  new  ideas  and,  for  this  reason,  the 
proof  is  given  in  some  detail.  Several  examples  are  given  Ln 
Section  3.  In  Section  A,  similar  results  are  presented  for 
equations  with  infinite  delays.  In  Section  3,  we  give  an 
example  in  two  dimensions  for  which  ‘■table  Mopf  bifurcation 
occurs  with  decreasing  delay.  In  Section  h,  we  give  an  in- 
troduction to  some  of  the  methods  available  for  nonautonomous 


equati  oris . 
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2.  The  Hopf  hi  furcation  theorem . One  of  the  simplest  vays 
in  which  nonconstant  periodic  solutions  of  autonomous  equa- 
tions can  arise  is  when  an  equilibrium  point  changes  from  being 
stable  to  being  unstable  as  parameters  vary  in  the  equation. 

A periodic  orbit  can  bifurcate  from  the  equilibrium,  the 
process  being  generally  referred  to  as  Ilopf  bifurcation.  For 
ordinary  differential  equations,  several  proofs  of  the 
existence  of  a Hopf  bifurcation  have  been  given.  In  some  way, 
they  all  involve  ultimately  the  implicit  function  theorem  and 
the  technicalities  of  the  proofs  are  minimal. 

In  the  statement  of  the  Hopf  bifurcation  theorem,  one  must 
always  impose  differentiability  conditions  on  the  vectoi  field 
in  order  to  obtain  smooth  bifurcation  curves.  For  ordinary 
differential  equations,  these  conditions  are  not  very  restrictive. 
For  functional,  differential  equations,  the  obvious  differenti- 
ability conditions  can  eliminate  the  discussion  of  variations 
in  important  parameters.  For  example,  consider  the  equation 

(2.1)  xft)  = A(6) x (t)  + B(6)x(t-r)  + C(6)x(t-s) 

+ f(6,x(t) , x ( t - r ) ,x(t-s)) 

where  e !R,  r > 0,  s > 0 are  considered  as  parameters, 

A(6)  . P (>')  , C ( 6 j are  n * n constant  matrices,  f(6,x,y,z)  as 
well  as  the  first  derivatives  with  respect  tc  x,y,7  vanish 
at  x = y = z = 0.  Suppose  the  characteristic  equation  of 


the  linear  part  of  Equation  (2.1) 
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(2.2) 


det  A(A,*,r,s)  = 0 

A ( A , <5 , r , s ) - At  A(6)  - R(6)e'Xr  - C(<5)e 


for  (A,r,s)  -•  ( 6 , s^)  has  a pair  of  purely  imaginary 

roots  jv(j,-iv  ,Vq  > 0,  and  ail  other  roots  have  negative  real 
parts . 

The  fundamental  prr  ’ - is  to  discuss  the  existence  and 

stability  of  small  nonconstant  periodic  solutions  of  (2.1) 
for  (6,r,s)  near  (6q,Tq,Sq)  and  which  vary  in  a smooth 
way  in  (<5,r,s,).  It  i.s  not  very  restrictive  to  assume  the 
right,  hand  side  of  Equation  (2.1)  is  continuous  and  has  a 
continuous  first  derivative  in  A.  One  can.  then  prove  that 
the  solution  of  Equation  (2.1)  is  also  continuously  differ- 
ent in  Mo  in  <5.  However,  regardless  of  the  space  of  initial 
data  for  Equation  (2.1),  the  solution  mapv.ill  not  be  differ- 
entiable in  (r,s).  At  first  sight,  this  makes  it  unclear 
how  to  solve  the  above  problem.  The  important  feature  that 
makes  the  problem  tractable  is  that  every  per  iodic  solution  of 
Equation  (2.1)  must  have  one  more  derivative  in  t than 
f(6,x,v,z)  has  in  x,y,z. 

The  Hopf  bifurcation  theorem  is  not  an  elementary  exercise 
for  functional  differential  equations.  The  verification  of 
the  hypotheses  necessary  to  apply  the  implicit  function 
theorem  uses  a number  of  special  identities  for  linear  systems 
with  constant  coefficients.  Since  we  are  going  to  state  the 


5- 


theorem  so  that  it  will  be  applicable  to  variations  in  the 
delays,  a proof  of  the  theorem  will  be  indicated  and  is  based 
on  the  proof  of  a less  general  result  in  [7,  p.  246]. 

Suppose  r > 0 is  a given  real  number,  U!  = (-''",'*•)> 
is  an  n- dimensional  linear  vector  space  over  the  reals  with 
norm  |*|.  C ( l a , b ] ^lRn)  is  the  Banach  space  of  continuous 
functions  mapping  the  interval  [n,h]  into  1Rn  with  the 
topology  of  uniform  convergence.  Tf  [a,b]  = [-r,0],  we  let 
C = C f ( - r , 0]  ,IRn)  and  designate  the  norm  of  4 by 
| <]>  | “ *>up  r<0<Q  I ^ ) I • Suppose  71  C 1R  is  an  open  set  (the 
parameter  space),  f : ,Q  x C -*■  lRn,  b:  1?  x C ->  !Rn  are  continuous, 
L(a)|>  j.s  linear  in  P , ff^,^)  has  continuous  first  and  second 
derivatives  in  <P , f(u,0)  = 0,  Df  (« , 0) /?<p  - n and  consider 
the  equation 


(2.?)  x(t)  - L (°0 x t * f(«,xt), 

where  x (9)  " x(t+9),  -r  < 0 < 0. 

Our  first  hypothesis  is  the  following: 

(lip  The  character istic  matrix, 

; (2.4)  A (« , A)  = A I - L(ot)eX'T, 

. 

where  f is  the  n * n identity  matrix,  .is  continuously 
diffei entiable  in  a,  there  is  a purely  imaginary  characteristic 
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root  ■ q ■ * v Q j vq  > 0,  for  a = «()  and  no  other  characteristic 
root  Vj  t X q . ^ g of  the  characteristic  equation 

(2.5)  det  A (ot , X)  --  0 


a a rt0  satisfies  X = mXp,  in  an  integer. 


Since  f(u,<£)  is  continuous  and  linear  in  f , Micro  is 
n y n matrix  function  qlct,1))  of  bounded  variation 
in  9,  -r  < 0 < 0,  such  that 


(2.6) 


h (°0 


0 

[dn(ct , e)  j q,  (o). 
-r 


Along  w ;i  th  the  linear  equation 


(2.7) 


x(t)  = L (a)x 


consider  the  formal  adjoint  equatio 


(2.8) 


on 


f0 

y(x)  = - y(x-0)dn(a,e) 

' -r 


and  the  bilinear  form 


0 r 9 

(2-9)  (4* ,<PD  = lH0)$(0)  - +(5-0)  rdn(a,0)]t(£)d$ 


- r'  0 


defined  Tor  't>  e C and  r e C f f U , r]  ),  JR“  space  of 
n-d  imens.i  onn 1 ro w vectors. 

lor  the  characteristic  values  (iv _ > - 5 v ) of  the  character 
istic  equation  det  £•(«  A)  = 0 corresponding  to  the  linear 
equation  (2.7)  for  a - a there  are  two  .l  inearly  independent 
solutions  h cos  v t,  b sin  v^t  foi  some  n -d  imens  i ona  I column 
vertm  b.  In  the  same  wav,  their  are  t'o  lineas.lv  independent 
solut  ions  c:  cos  vf)t,  c sin  v(.  x of  the  formal  ad  joint  equation 
(2.S)  for  a - a . In  Chapter  7 of  17].  the  following  remarks 
are  proved.  If 


<t>a^  = (T-,  ,<f>2)  » = b ccs  v0e 


= h sin  v 0 , -r  < 0 < 0 


<K(s)  = c cos  'xns 


0 lv. 


^2(S)  “ C sin  V’  0 < 5 r 


Of',,  T,1  Of',,  TO 

. * * . .L  X i.  t. 

t»a  »<f,a  ) * 

0 0 OfV'V 


where  ( , ) is  the  bilinear  form  in  (2.9),  then  t he  2 * 2 

* * _ i * 

matrix  (’!'  , 1>  ) is  nonsingular.  If  = (T  ,t  ) '7  , th 

o &o  0 0 u0 

(¥„  , ) = I and  one  can  decompose  C as 
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Thi  s 

decomposition  of 

C d e f i n e s z 

i project  ion 

r with 

0 

ion  is  made  in 

V 

b , , n-tu  )c 
n 0 

= Q • When 

n 

i decoinpos  i t 

thi  s 

way,  we  say  C 

is  decomposed 

by  the  set 

of  characteristic 

values  { i v , - i v . 


The  following  result  is  proved  in  exactly  the  same  way 
as  Lem  it!  a 2.2  :in  [7,  p.  17)]. 


Lemma  2_.  1.  If  (l!^)  is  satisfied,  then  there  is  a c > 0 and 
a simple  characteristic  root  X (a)  of  (2.5)  which  is  continuous 
together  with  its  first  derivative,  Ira  X («)  > 0,  J ct  - ct  { < 6, 

X = iv  . Furthermore , C can  be  decomposed  by  (A  (a) ,a  («) } 
as  C * t’a  fi)  Qa , dim  P = 2,  and  the  corresponding  projection 
operator  n((  is  continuous  together  with  its  first  derivative 
in  <* . 

With  this  lemma  and  f - fif.  ,<{>_)  a basis  for  P.  , it 

X .1  L Ct 

follows  that  there  is  a 2 * 2 matrix  B (a ) continuous  and 
continuously  differentiable  in  a such  that 


Ta(°)  = $a(0)exp  B ( a ) 0 , -r  < 0 < fl. 
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Also,  the  eigenvalues  of  B(a)  are  X (v)  ,!*(«)  . Furthermore, 
we  may  assume  by  a change  of  coordinates  that 


B (a)  vqB0  + B (a) 


(2.1.1  ) 


0 1 
1 0 


B,  =? 


(a-«pK  (•'« ) (!i-''n)  ■ V (a) 


- ('-t-«0)-Y  ('0  (u-a^  ) • (a)  J 


where  ?(•*), Y(ci)  e 1R  are  continuous  and  continuously  differ- 


entiable in  a for  |a-a^l  < 6. 


V.v  may  now  state  a generalization  of  the  Hopf  bifurcation  theon 


Theorem  2. 1 . Suppose  L(''0<j>,  f (u,*)  satisfy  (11^), 


(H?)  for  any  K > 0,  <P  e c with  d <t>/dd  r.  c,  |d'j>/dO|<  K, 


the  function  L(a)<j>  + f (a , <f ) has  a first  derivative  in  a 
which  is  continuous  in  a,$, 


(II.,)  C(a  ) i 0 


where  •;  (a)  is  defined  in  Equation  (2.11).  Then  there  is  an 

.1 


e > 0 such  that  for  a e F , |a|  < e,  there  is  a C -manifold 
ra  el?  of  cod iiiiens ion  1,  I’a  is  continuous  and  continuously 


differential  in  a 


1’  = (a  tf  : Re  X(a)  = 0,  |a-M  | < e} 


where  X(ot)  is  given  in  Lemma  2.1,  such  that  for  every 


J 


a t r there  ;is  a function  to  (a , a ) . an  <o(»,a)  -periodic 


function  x"  (•<  ,a)  continuous  tpgcthe r_wi  th  thcii  _ first. 


derivatives  in  a,  a,  ui(a  , 0)  = w - 2«/v  , x fa,0) 


0, 


xq  (!a , a)  * 't>u  col  (a  , n)  + of  | a 


as 


0 


and  x'f't,a)  is  a solution  of  Equation  f2.31.  1 ur t hc_t more . 


for 


a i - w , 


< c every  id  - n er  iodic  solution  of 


Equation  (2.3)  with  |x 


< c 


must  bey  o f the  above  type 


except  for  a translation  in  phase, 


Before  proving  the  theorem,  let  us  discuss  some  implications 
If  k - 1,  that  is,  a o 111,  and  A (a)  - |j(e)  + io(u)  is 
tlie  characteristic  root  given  by  Lemma  2.1,  then  (II,)  is 
equivalent  to 


d>(u{)) 


Hci  “ 


t 0. 


Thus,  Theorem  2.1  is  the  usual  Hopf  bifurcation  theorem.  In 
fact,  for  k = 1,  the  set  can  be  described  by  a smooth 

curve  u(o).  Thus,  for  jaj  < c.  , there  is  a function  a(a) 
continuous  together  with  its  first  derivative  such  that  a(0)  = 
and  x'  (a)  = x*(«(a),a)  is  an  to  (a ) = w (ot  (a ) . a)  periodic 
solution  of  Equation  (2.3).  The  bifurcation  curve  oof  a)  is 
shown  schematically  in  Figure  1. 

If  k - 2,  then  the  bifurcation  surface  in  the  (a, a) -space 
is  shown  schematically  in  Figure  2.  For  cadi  a,  there  is  a 


Figure  ] Figure  2 

smooth  curve  of  a values  for  which  there  exist  nonconstant 
periodic  solutions  of  Equation  (2. .I).  Tf  k > 2,  the 
picture  is  similar  except  F is  a surface  of  codimension  ] 

id  k 
i n il?  . 


To  prove  Theorem  2.1,  we  first,  obtain  the  bifurcation 
equations  for  the  small  amplitude  periodic  solutions  of  the 


- 12- 


equation  (2.3)  foi'  |a-a  | < 6 which  have  n period  close 


to  ?.v / \'  As  in  ordinary  differential  equations,  we  normalize 
the  period  to  2?r.  Let  B e (-1,1),  = ?n/v  t " (1-*B)t, 

x(t‘-(l ) ■ u(T  + O/(l+0))  , -r  < (3  < 0 and  define  ur  as  an 

1 7 1 

clement  of  the  space  C([-r,0]^TRn)  eiven  by  ur  ( - ) - 

1 > p 

u ( t *• .' / ( ;|  + p,) ) , -r  < B < 0.  Equation  (2.3)  is  then  equivalent 


(2.12) 


" L- = Cl- r>)  [h(«)uTfS  + na,u.tjf5)) 


If  this  equation  has  an  -periodic  solution,  then  'Equation 
(2.3)  has  a (1  + 6) - per  iod  ic  solution,  and  conversely. 
Equation  (2.12)  may  lie  rewritten  as 


(2.13) 


on11 ' u«0)u,  * 

N(6,a,uT,utj6)  > - I.(«0)uT  ► (1*  6)  f (« ,u  t _ p) 


We  consider  Equation  (2.13)  as  a perturbation  of  the  autonomous 
linear  equation 


(2.141 


- r.<V“. 


We  know  that  the  columns  of 


U(t)  ur  (P)  exp  R(ct  )t,  x 
0 u 


* 
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form  a basis  for  the  -periodic  solutions  of  liquation  (2.14) 
and  tbo  rows  o 1' 


V(i)  d = f[^p  - B(an)-t)Y  (0),  t e JP. 

0 


form  a basis  for  t/e  u>g- peri  odic  solutions  of  the  formal 
adjoint  liquation/  (2.8). 

need  the  following  lemma  from  [7,  p.  209]. 


Lemma  2.2,  Let  9*^  be  the  Banach  space  of  cont i nuous 
oig -periodic  functions  with  values  in  TRn  and  the  topology 
of  uniform  convergence. lor  any  I e , rhe  equation 


(2.15) 


“ItT1  = L<«0)ur  ' f m 


has  a solution  in  if  and  only  if 


(2.16) 


u) 

r 0 


y ( t) f ( t) dt  - o 


for  all  flp -per iodic  solutions  of  the  formal  adjoint  equation 

(2.8).  Furthermore,  there  is  a projection  J:  9'  -*  9* 

W0  ^0 

such  that  the  set  of  f in_  9*^  satisfying  (2.16)  i_s 
( I - J ) 9>^  anti  there  is  a continuous  linear  operator 
K:  ( I - J)  9ifa  -►  ( I - tv)  9>^  such  that  Kf  is  a solution  of 
(2.15)  lor  every  f e (I-.J)  9>^  . The  operator  rr  is  any 
projection  onto  the  a -periodic  solutions  of  (2.14). 
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1 he  opera  Lor  J is  given  by 


(2.17) 


.]£  = V'  [ V (s) V ' (s) cts] 

1 0 


V ( s ) f ( s ) d s 


where  V'  is  the  transpose  of  V. 

This  lemma  simply  states  necessary  aiul  sufficient  con- 
dition? for  the  existence  of  periodic  solutions  of  non- 
homogeneous  linear  equations.  Knowing  these  conditions,  one 

can  always  subtract  off  part  of  any  f e & , namely  Jf,  so 

0 

that  1 he  orthogonality  conditions  (2.16)  are  satisfied  for 
(I-.I)f.  Thus,  a solution  will  exist.  We  can  always  choose 
a particular  solution  to  he  orthogonal  to  the  solutions  cf 
the  homogeneous  equation.  This  is  K(I-J)f. 

Using  Lemma  2.2,  we  Know  that  every  <»>  -periodic  soLution 
of  Equation  (2.13)  except  for  a translation  in  phase  is  a 
solution  of  the  equations 


(2.18a) 


(2.1810 


u(-)  = U ( • ) col (a , o)  + K(I-J)N(6,u,u< ,u>  g) 


JN ( ft , « , u ,u  J = 0 

• * 9 F 


and  conversely. 

One  can  now  apply  the  implicit  function  theorem  to  solve 
Equation  (2. 18a)  for  u = u'v(a,ft,a)  for  a, ft, a in  a 
sufficiently  small  neighborhood  of  zero,  u*(a,0,0)  - U(Ocol(a,0) 
= o(|a|)  as  [a|  -*■  0.  The  function  u is  continuous  together 
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with  derivatives  up  through  order  two  from  an  application 
of  thr‘  implicit  function  tlieorem.  However,  without  some 
additional  information,  we  cannot  prove  it  is  differentiable  in 
ft, a.  Since  u*  (a,  (t)  satisfies  (2.38a),  it  is  aut.o- 

mat.icn.lly  continuously  differentiable  in  t with  hounded 
derivative  for  |a-ct^|  < <5.  If  we  use  this  fact  and  the  form 
of  the  implicit  function  theorem  given  in  Lemma  2.7  in  (7,  p.  2361, 
one  obtains  the  following  important  result. 

Lemma  2.3.  If  (11^)  , (lip  are  satisfied,  then  the  function 
u * ( a , f , a ) j_s_  conti nuous  t oget he r with _i ts  f i rst  d er  ivatives 
in  a,S,u  and  second  derivatives  in  a. 

from  the  above  lemma,  all  m -periodic  solutions  of 
Equation  (2.13)  are  obtained  by  finding  the  solutions  a,B,ct 
of  the  bifurcation  equations 


(2.19) 


JN(B,«,u"  (a,  3, a)  ,u*  (a,  B,cx)) 

• • y P 


0. 


Using  the  definition  of  J in  Equation  (2.17) 
(2.19)  is  equivalent  to  the  equation 


w 


def 


(2.20)  G(a,B,c<) 


0 -B(o  )s 

o V II 


0 


l'a  (°)N(B,a  ,u"(a,fi,m),us^g 


the  Equation 


(a,B,a))ds  = 0 


Prom  the  above  discussion,  it  follows  that  it  remains  to  solve  the 
equation  G(a,8,«)  = 0,  which  represents  two  equations  in  the 
parameters  a,  (3, a. 


L 


Since  G(Oy|?,*)  = 0 for  all  P,a,  define 


(2.7-1)  II(a,P,a)  = G(n,P,«)/a. 

This  function  is  continuous  together  with  its  first  derivatives 
in  a , (5 , a . 

The  nonzero  small  amplitude  periodic  solutions  of 
Equation  (2.3)  are  given  as  x(t)  - u*  (a  , 6 ,;0  ( t/ 1 + P)  ) where 
a,P,«  satisf>  the  equation 

(2.22)  II(a,3,«)  = 0. 


Our  next  objective  is  to  compute  the  Jacobian  of 
H(a,6,o)  with  respect  to  6,u  evaluated  at  (a,3,»)  = 

(0,n,a  ).  From  the  definition  of  G in  (2.20),  N in  (2.13) 
and  u'  one  easily  observes  that,  for  e^  - col(l,0), 


H ( 0 , 0 , ct ) 


/J0  -B(a  )S 

e T (0){l,(a)U  e - I,(P)U  e.}ds 

- 0 0 5 L s 1 


( 0 -B(a  )s 

e fa  (0)(L(a) 

- o n 


b(%)s 

h(an)  ) J,.  e ' e.ds 
0 L 


since  U - $ exp  B (a _) s . 


We  need  tlie  following  Lemma  3.9  in  (7,  p.  179]. 
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Lenima  2.4.  Suppose  the  conditions  o f Lemnu  2.1  arc  satisfied. 


Then 


dR(a  ) dL(<0 

(L  . -s  (o')  -_-u~  • 

du  % d*.  0 


Using  Lemma  2.4, 


ail(0,n,an)  f“l)  -M«0)s  clB(a  ) U(a  )s  dh(u  ) 

0 V - e ° — zr0~  e *jds  = — g*  — e1 

n 


(2.2.3) 


',-n»0) 


where  C(<’0)  , v(k„)  are  defined  in  (2.11). 

(0  compute  tl\e  derivative  with  respect  to  1?,  we  note 


first  that 


H<°*P»V  " J0  e 


0 - B (a q)  s 


(°){(1+B)L(ao)Us,Bel  ' L(VUsel)ds 


Writing  this  as  two  separate  integrals,  changing  s into 
s/(l+P)  ’n  the  first,  integral  and  noting  that 


dU(s/fl+B))  _ T rn  -)  II 
“T(s7lT+ltj  L 1 V US/  (1  * B) , e 


one  sees  that 


0 -B(a0)s  B fo.  ) «: 

H(0,Bf‘O  - p e f (0)  * (0)e  0 B(an)e.ds. 

J0  ao  0 0 1 

If  x is  .i  solution  of  liquation  (2.14)  and  y is  a solution  of 
the  adjoint  equation,  then  (yt,x  ) = constant  for  all  t 
There fore . 


B(a  )(>•) 

(e  1 ro  (o),  t f 0)o  u ) 

0 0 

*B(a  )s  BfuJs 

e J T (0)t  (0)e  " 

0 0 


f o -B(an)(s+e-0)  B(a  )(S  + C) 

e U I'  (O)d.i(a  0)t  ( 0 ) e 1 

J-rJn  o o 


for  nil  s e 1R.  Integrating  this  from  0 to  w and  using 
the  fact  that  the  second  integral  is  zero,  one  obtains 


H(0,6,«0)  = B"0B(a0)e]  = -2*3  col (0,1) 

(2.21) 


3H(a,B,a0)  f 0 


If  H = col  (IL  , H,)  then  Relation  (2.21)  implies 


311,(0, 0, a ) 


il  (0,0, aQ)  = 0,  -~-W 


Conscquently , the  implicit:  function  theorem  implies  there  is  an 


e > 0 and  a function  B (o,a)  continuous  together  with  its 
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first.  derivatives  in  a,«M  , | a-ct^  | < c,  ( a | < c,  P*(Oq,0)  = 0, 
such  that  P#'(w,a)  satisfies  the  equation 

(2.2  9)  !!2(a,B*(ala)fu)  = 0 

and  is  the  only  solution  for  |(?|  < c,  |a-u  | < c,  |g|  < e. 

Up  to  this  point,  we  have  the  following  lemma. 


Lemma  2.5.  If  (llj ) , (It-,)  are  satisf  ied,  then  there  is  an 
£ > n»  fillictions  3*  (a, a),  u*  [a , B,a)  , continuous  together  with 
their  first  derivatives,  u*  (a,  £>,«)  to  -periodic,  |aj  < e. 

I a ■ u q I < e,  6 (up»0)  = 0 , n " ( 0 , 0 , ot  ) = 0,  such  that  for 
|“-")0  I < e,  | cv  ~ ct  0 I < er)  > Equation  (2.3)  has  an  oi- per  .iodic 
solution  x with  |x|  < e if  and  only  if  x = u* (a , 6* (a , a) , a) 
except  for  a translation  in  phase  where  a, a satisfy  the 
bifurcation  equation 

(2.2.6)  h(a,a)diCH1(a,  p*(ct,a)  ,a)  = 0. 

.t. 

The  function  u (a,  6, a)  satisfies  (2.18a)  and_  fl  (a, a) 
sat  i s i i ejs  (2.5)  . 

Remark  2.1.  It  is  clear  from  the  proof  of  the  above  results 


that  if  we  assume  I,(u)$  ^ f(a,<f)  lias  k derivatives  with 

respect  to  f which  are  continuous  and  Hypothesis  (il7)  is 

t"  h 

satisfied  for  the  k'  derivatives  with  respect  to  ct 


when  '}•  has  k continuous  derivatives,  then  the  function 
h(u,n)  has  continuous  derivatives  in  u , a up  through  order  k. 

From  the  expressions  for  the  partia]  derivatives  of  II  with 
respect  to  «.£  at  (0,0,a^),  one  easily  determines  that 

h(«,a)  -■»  S(|(n)  + (u-an)^1(n)  + o(|u-a0|) 

(2’?''  ^(0)  = 

where  £(ut„)  is  given  in  (2.11). 

t * 

Proof  of  Theorem  2.1.  If  r = f«:  h(<*,a)  “ 0,  | a j < e}, 
then  the  Implicit  Function  Theorem  implies  the  assertions  about 
T if  It  - l . For  k > I,  it  .is  the  transverse  l i ty  theorem 
(see  (I,  p.  45]).  Then  u)(a,a)  = (1  + 3 ’('*>,a))w  , x"(«,a)  = 

£ :*r 

u' (a.p  fu,a),«)  satisfy  all  of  the  properties  stated  in  the 
theorem.  This  proves  Theorem  2.1. 

Theorem  2.1  uses  only  knowledge  of  the  linear  operator 
L(a)  in  liquation  (2.3).  We  know  nothing  about  the  specific 
structure  of  the  sets  T . To  obtain  a more  complete  picture 
of  the  bifurcation,  one  must  consider  nonlinear  terms.  bet 
us  suppose  bf'04*,  f («,<)>)  satisfy  all  conditions  above  with 
the  differentiability  conditions  up  through  order  k > 3.  Then 
the  function  h(«,a)  in  Relation  (2.27)  satisfies 

h(a,a)  = h (a , 0)  + p(<*)a2  + o(|a| '), 
h(a,0)  = (a-a0)  + C(a.Q)  f o(|a-a0|), 


(2.28) 


as  I a | + o , a -y  a ^ . Let  ns  assume  that: 


(2.29) 


■<V  * " 


The  function  h(«,a)/p(a)  will  then  have  a minimum  at  a = 0 
given  by  h(o  , 0) /p(&)  . The  bifurcation  equation  will  have  a 
solution  if  and  only  if  this  minimum  is  1e°-s  than  or  equal  to 
zero  and  will  have  no  solution  if  this  minimum  is  0.  The 
minimum  equal  to  zero  is  therefore  a bifurcation  surface  in 
the  a- space.  This  is  given  by 


(2.3(0  0 - M“jfl  . p-J-T  < oc|«.«0i) 


as  a <■  os 

0 

Under  hypotheses  ( 1 ) t liquation  (2.30)  defines  a smooth 
hyperplane  Q in  a neighborhood  of  a in  .1R  . Jn  a 
sufficiently  small  neighborhood  of  the  point  (a  ,0)  in 

x C,  on  one  side  of  this  hyper plane  there  are  no  nonconstant 
periodic  solutions  and  on  the  other  side  there  is  a unique 
noncoiu  I an t periodic  solution.  For  k ~ 2,  the  bifurcation 
diagram  is  similar  to  the  one  shown  in  Figure  3.  The  curve 
is  the  intersection  of  the  surface  in  Figure  3 with  the  plane 
a = 0 . 

The  direction  of  bifurcation  is  determined  by  p(&g). 

These  results  are  briefly  summarized  in  the  following 
result. 


addition , ( 1 1 ^ ) is  satisfied  and  p ( a ) f-  0 where  p ( a ) is 
given  in  (2.28),  then  there  are  neighborhoods  V Cli*1  o_f  ct^ , 

U in  C of  zero  and  a smooth  hyperplanc  PCI!  such  that 
V\T  = A1  U \2  and 

(i)  a e Aj  implies  no  nonconstant  periodic  solution  in  U 
(ii)  a e A2  implies  a unique  nonconstant  periodic  solution 


in  U . 
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A 1 s s general  version  of  the  lix  i r-  tencn  Then  v cm.  2.2  was 

obtained  by  Ru  i z-Cl  aeysson  (22,  Theorem  ft.''],  lie  considered 

1 

the  equation  (2.3)  in  the  space  IV  ’ in  order  to  obtain  the 
cent'  r manifold  and  then  applied  techniques  similar  to  the 
ones  employed  by  Chafee  f 4 ] . In  addition,  ho  was  able  to 
extend  the  analysis  in  Hale  [9],  llausrath  1111  to  obtain  the 


stab  it  ity  (if  t he  bifurc.nl  ion  when  ( I F ^ ) was  i\rcii"  t honed  in  say 
that  alj  characteristic:  values  have  negat  ive  real  parts  except 
the  two  on  the  imaginary  axis.  The  stability  properties  are 
determined  by  the  change  in  stability  of  the  zero  solution  as 
is  considered.  One  can  actually  improve  the  results  in  [22] 
as  well  as  extend  them  to  the  case  when  liquation  (2.3)  is 
considered  in  C.  More  specifically,  one  can  prove  the  follow- 
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3.  Examples  o.l  variations  in  the  delays.  Consider  the  .linear 
scalar  equation 

(3.1)  x = - jx(t-r^)  j x ( t - r ? ) 

where  r.j,r7  are  positive  constants.  The  specific  co- 
efficients :i ii  liquation  (3.1)  are  chosen  in  order  tc  male  the 
computations  simpler.  'Ihe  characteristic  equation  for 
Equa  t i on  (3.1)  is 

1 'Xrl  1 ~Xr? 

(3.2)  A + i e 1 ->  j e z = 0. 

To  dismiss  the  Hopf  bifurcation  for  a nonlinear  perturbation 
of  Equation  (3.1),  one  must  determine  the  curves  in  (r^ji^)- 
space  at  which  the  roots  of  Equation  3.2  have  zero  real  parts. 
Ruiz -Claeyssen  [22]  has  analyzed  this  problem  completely. 

There  are  infinitely  many  curves  for  which  there  are  purely 
imaginary  roots,  but:  we  are  going  to  concentrate  on  the  curve 
which  lias  the  property  .if  (r^.r^)  r.  T,  then  there  is  a pair 
of  purely  imaginary  roots  and  all  other  roots  of  Equation  (3.2) 
have  negative  real  parts. 

For  - r,  = r,  Equation  (3.2)  is  the  familiar  equation 

(3.3)  A + e'Xr  « 0 

There  is  a unique  r^  = ir/2  such  that  Equation  (3.3)  lias  all 
roots  with  negative  real  parts  for  0 < r < and  a unique 


v 
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pair  of  purely  imaginary  roots  1 ~ i i lor  r ~ i . The 
curve  I'  that  interests  us  is  the  one  in  the  ( r . , t ,) -plane 
which  passes  through  the  point.  (1|/2,'n/?). 

If  A - j.y  in  liquation  (3.2),  then  ve  must  have 

Y(r  +r  ) Y(r  -r  ) 

0 - cos - cos  

( 3 . '1 ) . , 

'IVY  Tlll'r2) 

Y = sin  ~2 cos — 

Since  y - 0 is  not  a solution  of  Equation  (5.2),  it  follows 
that  cos  (r,-r_)/2  f 0.  Thus,  we  must  have  y - tr/(r  + r ), 


T! 


COS 


’(yy) 

-*r1+rT 


Figure  1 

The  curve  F is,  therefore,  the  solution  of  this  equation 
and  is  shown  schematically  in  Figure  4. 

Now  consider  a perturbation  of  Equation  (5.1) 

x(t)  = - | x(t-rj)  - | x(t-r2)  + f(x(t),x(t-r1),x(t-r2)) 
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where  f(x,y,z)  is  continuous  together  with  derivatives  up 
through  order  three,  f ami  its  first  partial  derivatives 
vanish  at  x r'  y = z = 0.  Theorem  2.1  implies  there  is  a 
I lop  f bifurcation  at.  ever)’  point  in  f. 

To  obtain  more  information,  one  must  assume  something 
about  the  nonlinearity  f.  In  particular,  consider  the 
equn  t i on 


(3.5)  x(t)  = - y x(t-rj)  - j x(t-r?)  + Px’(t),  B f 0. 

If  wc  let  ctfi  = (r^.r^)  f-  T,  then  the  characteristic  roots 

X (o  ) , \ (ot n ) on  the  imaginary  axis  are  given  by  i vn  , - i)  q , 

Yq  - 1r  / ( r ^ + r ^ ) . Choose  bases  1(^1  ,T  (u())  ns  in  the  previous 
section  for  the  decomposition  of  C as  C ~ P(J  61  Qa  foi 
the  characteristic  roots  { i.y() , - i Iq  } . In  [2  2],  it  is  shown 
that  pfa^)  in  (2.28)  is  given  by 

(3.6)  p(af))  " '5cc(ap)fi 


where  k is  a positive  constant, 


(3.7) 


c(«0)  = (1-Y0  Sin  YqI-J/?-) 


i r 

v c k 2 


c (a  0 ) i 0, 


Theorem  2.2  implies  there  is  a neighborhood 
such  that  V\F  = U A,  and 


of  u e r 


-27- 

(i)  tlierc  is  no  small  nonconstant  periodic  solution  of 
(3.1)  for  (i j , r?)  c A j 

(ii)  there  is  a unique  small  nonconstant  periodic  solution 
of  (3.1)  for  (r1,r2)  e A, 

that  is,  a small  segment  of  r is  a bifurcation  curve. 

Rui.z -Cl  oeyssen  [22]  also  showed  that  c(fin)3  > n implies 
the  periodic  crb.it  is  asymptotically  stable  and  the  bifurcation 
occut s as  T is  crossed  from  right  to  Jeft.  If  c(u^)6  < 0, 
the  periodic  orbit  is  unstable  and  the  bifurcation  occurs  as  1 
is  crossed  from  left  to  right.  There  are  some  points  on  T 

where  c(c'q)  ~ 0.  No  information  is  available  near  these 
points  since  higher  order  approximations  in  the  bifurcation 
equations  will  be  needed.  However,  we  do  sec  that  a two- 
parameter  problem  is  very  different  from  a one  parameter  problem. 
The  direction  of  bifurcation  can  change  as  we  move  along  r . 

Stech  [2d]  lias  also  observed  this  phenomena  in  problems  with 
two  delays. 

An  obvious  generalization  of  Equation  (3.1)  is 
(3.R)  x(t)  ~ -(l-«)x(t-rj)  - ux(t-i^),  0 < u < 1, 


whose  characteristic  equation  is 


-Xr 


1 


(3.P) 


X + (1  -u) e 


* ae 


0. 
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The  analysis  of  Equation  13.9)  lor  all  « r (0,1],  i'  > 0,  > 9 

is  difficult.  From  our  point  of  view,  we  would  like  to  know  for  a 
fixed  ci  the  curve  V in  (r^  , r ^)  - space  which  passes  through  the 
point  (n/2,n/2)  and  lias  the  property  that  (r^.r^j  e r implies 
two  roots  purely  imaginary  and  the  remaining  roots  with 
negative  real  parts.  This  curve  changes  in  a complicated  way  as 
a function  of  a and  is  not  completely  understood. 

Hie  equation  for  the  purely  imaginary  Toots  \ = iY  of 
Equation  (3.9)  are 

(a)  (l-a)cos  YTj  + a cos  yr_,  = 0 
(3.10) 

(hi  (1 -a)  sin  yr.  + a sin  yr„  = y. 

1 i. 

For  a = o,  the  solutions  of  this  equation  are 

yrL  = (2k+l) n/ 2 , y = (-l)k,  k - 0,1,2,..., 

and  the  special  curve  T is  given  by  r^  ~ ti  / 2 , r^  arbitrary, 
and  corresponds  to  k = 0,  y = 1. 

For  a = 1/2,  Equation  (3.10)  can  bo  written  as 
Equation  (3.4)  and  has  the  solutions 


= (2k+l)  7i/  (r1  + r2) 


(3.11) 


(2k+l  )7i 


= cos  (2k+lW 


rl+r2 


r2~rl 

rdf, 


k = 0,1,2, 


As  remarked  earlier,  the  special  curve  F corresponds  to  k = 0 
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and  is  shown  schematically  in  Figure  '1.  To  understand  how  F 
changes  with  c* , it  is  necessary  to  investigate  the  solutions 
of  liquation  (3.11)  for  other  values  of  k.  for  k = 1;  that  is, 
y = 3t/ ( r 1 + r 2 ) , 


T 1 


(3.12) 


3 n 3 1 2 L 

— - — cos  — — - 

rl  r2  2 r2  1 


The  solutions  of  liquation  (3.12)  are  depicted  in  Figure  5a  for 
r~  > t-j  . in  this  figure,  we  have  also  included  the  original 
curve  F . 

The  value  of  y = (2  k + 1 ) fi/  (r  ^ + r ^ ) , k = 0,1,  approaciics  zero 
as  r.,1",  approach  infinity  along  either  cf  these  curves; 
that  is,  for  large  values  of  r.  + r ^ , there  are  many  purely 
imaginary  roots  of  liquation  (3.11)  near  zero.  Consequently, 
for  a small  change  in  a from  a = 1/2  in  liquation  (3.10),  it 
is  conceivable  that  the  curves  representing  the  solutions  of 
Equation  (3.10)  for  large  values  of  r.  + r.,  can  have  a 
different  topological  structure  than  for  a « 1/2.  This 
suspicion  was  verified  on  a computer  by  M.  Michaud  and  T.  Lewis 
for  the  values  of  a depicted  in  Figure  5 and  corresponding  to 
the  two  curves  associated  with  k = 0,1  for  a = 1/2. 

In  Figure  5b,  the  smooth  parabolic- 1 ike  curve  corresponds  to 
the  curve  k = 1 , y - 3 77/ (r^M^)  for  a = 1/2  and  the  other 
curve  to  k - 0 , y = ^/(r^+r^)  for  a = 1/2.  At  the  cusp 
point,  there  are  two  distinct  pairs  of  purely  imaginary  roots  on 
the  imaginary  axis.  The  number  of  roots  with  positive  real  parts 
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near  this  cusp  is  indicated  in  Figure  5b.  The  llopf  bifurcation 
theorem  can  be  applied  along  each  smooth  branch  passing  through 
this  cusp.  The  end  result  is  that  it  is  possible  to  have  the 
same  number  of  periodic  orbits  as  roots  with  positive  real  parts. 

For  r.,  = 3r. , it  is  not  difficult  to  show  (see  Stech  [24]) 
that  :in  a neighborhood  of  a = 1/4,  liquation  (3.10)  lias  one 
solution  for  u < 1/4  and  three  solutions  for  ct  > 1/4  with 
yr^  near  n/2.  This  is  a different  type  of  phenomena  titan  the 
one  exhibited  near  the  cusp  where  the  purely  imaginary  parts 
were  different  depending  on  which  branch  was  followed.  For  a 
near  1/4,  the  curve  is  smooth  but.  lias  a "kink"  in  it  (see 
Figure  5d)  . Near  o.  - 1/4,  = 3r^  , it  is  possible  to  obtain 

three  periodic  orbits.  Similar  phenomena  have  been  observed 
by  Nnssbaum  121). 

For  larger  values  of  r,  than  shown  in  Figure  5,  more 
"kinks"  and  cusps  must  appear  because  of  the  fact  that 
Equation  (5.10)  are  almost  periodic  functions  in  r^,^. 
Therefore,  similar  bifurcation  phenomena  occur. 

These  results  are  important  because  they  have  the  follow- 
ing implications.  If  the  ratio  of  the  delays  is  constant,  say 
r 2 = flfj , 0 > 1,  then  it  is  possible  to  have  an  a and 
numbers  1 < 0n  < 6^  such  that  along  the  line  - Pr^ 

(i)  jf  1 < 0 < pg  there  is  a unique  intersection  with 

V 

(ii)  if  |3q  < 3 < , there  are  several  intersections 

with  r0, 
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Mi  i. ) if  3^  < ft,  there  is  a unique  intersect i on  with 

V 

Thus,  there  may  be  more  than  one  bifurcation  of  a nonlinear  | 

perturbation  of  Equation  ( 3.8)  as  we  move  along  a ray 

f2  ■ Bi  t - j 

Tor  the  particular  problem 

i 

t 

I 

(3.13)  ::(t)  = - ( (1  -a)x(t.-r1)  *•  u>:  (t  - r ,) ) [ Mx  (t)  ] 

0 < ex  <r  J,  Nussbaum  [21]  lias  obtained  some  very  interesting 
global  results  for  r]>r?  not  necessarily  near  the  curve  ra 
but  being  in  a wedge  - Br^  , 1 < B < 3.  The  technique  is 

to  use  ojertivo  fixed  point  theorems  and  he  has  proved  the 
following  results: 

(i)  if  1 < p < 2,  r?  ~ gr^  and  rf  is  sufficiently 
large,  there  is  a nonconstant  periodic  solution  of  Equation 
(3.11)  . 

fii)  if  2 < p < 3,  r,  - Br. , there  aie  some  ranges  of  u 
such  that  for  all  r ^ sufficiently  large  there  are  at  least 
trvo  solutions  of  Equation  (3.13). 

Our  results  on  local  bifurcation  theory  Avjll  apply  equally 


as  \i c 1 1 to  the  equation 
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One  cnn  study  the  bifurcation  as  a function  of  , r . One  of 
the  nni.li  difficulties  always  consists  of  analyzing  the 
c ha  r a c t e r i.  s tic  e qua  t i o n 


X 


+ a 


’ " l in 

k(0)e  dO 

-1 


0. 


The  existence  of  a bifui cation  can  he  then  determined.  The 
specific  nature  of  the  bifurcation  depends  on  the  nonlinear 
terms.  As  remarked  earlier,  the  most  efficient  scheme  for 
treating  the  nonlinear  terms  near  x - 0 is  the  method  of 
averaging  (see  [5]). 

Some  global  results  for  Equation  (3.14)  have  been  re- 
cently obtained  by  Walther  [27]  for  |r-l|  small. 
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4.  Infinite  delays.  In  this  section,  we  summarize  part  c£ 
the  Miosis  of  Lima  [161  on  the  I lop  I'  bifurcation  theorem  with 
infinite  delays.  In  the  case  of  finite  delays,  the  solution 
x of  a retarded  equation  will  belong  to  ('  for  t > r 
regardless  of  the  space  of  initial  data  for  which  one  defines 
solutions.  This  is  the  reason  why  Hie  initial  space  C cr 
JRn  x l.P(  { - r , 0]  ,IRn)  or  V-' 1 ’ p ( [ - 1 , 0 ] ,IRn)  has  little  effect 
on  the  qualitative  theory.  For  infinite  delays,  the  space 
of  initial  functions  becomes  more  crucial  since  all  of  the 
past  history  is  contained  in  x . In  the  first  part  of  this 
section,  we  consider  a special  case  of  the  initial  space  in 
Lima  110]  and  later  make  some  remarks  about  other  spaces. 

Suppose  g,  G:  (-*>,0]  -*■  JR,  are  nonnegative  continuous 
functions  satisfying 


(4.1) 


g(t+s)  < G(t)g(s) 
f0 

g < °° 

-oo 

there  is  a e (-"°,0]  sucli  that  G(Bq)  < 1. 


For  example,  if  g(0)  = exp  A0,  then  G(0)  = exp  AO  and  all 

- 2 

hypotheses  are  satisfied  if  A > 0.  The  function  g (0 ) = 0 
will  not  satisfy  the  hypotheses.  Let 
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(4.2)  I’^  “ {equivalence  classes  of  all  measurable  functions 

'!> : f -co , 0 J -►  ]Rn  such  that.  1 1>  | f | ^ ( 0 ) j ^ + 

+ f EC0)  |4»(3)  | 2 cl  0 } 1 7 2 < co}. 

~ cn 

Suppose  0 C 3R^ , f , I.:  ft  x B7  ->  lRn  are  continuous,  1-  Ca ) ‘J* 
is  linear  in  4',  f(u,<{>)  has  continuous  derivatives  in  4 up 
through  order  two,  £(«,(')  = 0,  Df(u,0)/;i4  ~ () . Consider  the 

lineal  equation 


(4.3) 


x(t)  = L (a ) x t 


as  well  as  the  perturbed  equation 


(4.4) 


x(t)  = L (ct  \ x t + f(a,xt) 


Since  L(a)<j>  is  linear  in  4,  it  can  be  written  as 


(4.5) 


L(a)<j>  = A(a)4(0)  + k(a  ,9)4 (6)g(6)d9 


where  k(a,*)  t [ ( -<*> , 0)  ,TRn  ] and  A (a)  is  an  n * n 

constant  matrix.  In  the  remainder  of  the  section,  we  assume 
L(a)j>  is  given  by  Relation  (4.5)  and  that  A(u),  k(a,‘)  are 
continuous  in  a.  The  characteristic  equation  for  Equation 
(4.3)  is 
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(4  . 6) 


dot  A fct , A ) ==  0 
■ 0 


A(cx,\)  A I - A (ci) 


k(a,0)eX,ifi(e)dO 


The  formal  adjoint  equation  is 

(4.71  y(t)  - -y(O)Afa)  ■ 


y(t-O)k(«,0)c(0)dO 


and  the  corresponding  bilinear  form  is 

■0 


(4.p)  - nomo)  + 


0 


c/-' ( 0— u )k('-i  ,0  )i(u)du 


g f 0 ) d0 


defined  for  9 e B.?  and  an  n-diraensional  row  vector  function 
in  some  appropriate  class  of  functions.  The  specific  class  is 
not  important  to  us  hero  and  may  he  found  in  [16].  It  is  only 
necessary  to  remark  that  <p  ,<j>>  will  be  meaningful  for  all 
functions  considered  below. 

Before  stating  the  Hopf  bifurcation  theorem,  let  us  point 
out.  the  important  implications  of  G(Bq)  < * n Relation  (4.1) 
If  x(t)  is  the  solution  of  Equation  (4.31  for  <f>  e B^  and 
T ( t ) j x (<}>),  t > 0,  tlien  (T(t),  t > 0}  is  a strongly 

continuous  semigroup  of  linear  operators  on  B.,.  If 
G(Pn)  < 1,  then  there  is  a y > 0 such  that  the  essential 
spectrum  of  T(t)  belongs  to  the  circle  of  radius  e u , 
t > 0 (see  [16],  [19]  or  [10]).  Therefore,  only  elements  of 
the  point  spectrum  of  T(t)  lie  outside  this  circle  and  these 
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the  ('harac.t  cr  Lit.  ic  Equation  (4.6).  This  makes  the  llopf 
bifurcation  theorem  feasible  since  one  can  have  two  character- 
istic values  on  the  imaginary  axis  ami  know  that  it;  is  possible 
to  have  the  other  elements  of  the  spectrum  of  T(l)  inside 
the  uni t circle. 

I'll?  discussion  below  is  a slight,  genera  1 i za  t i on  of 
Lima  [ I6|  taking  into  account  the  fact  that  we  know  from 
Section  2 that  one  can  prove  bifurcation  theorems  vitheut 
having  L (ci ) d,  f («  ,<P  ) differentiable  in  i for  all  functions 
if  we  are  very  careful  in  the  analysis. 

The  hypotheses  imposed  for  a bifurcation  theorem  are  the 
same  as  the  ones  in  Section  2 for  finite  delay. 

(IIj)  The  characteristic  matrix  A (a  ) in  Relation  (4.6) 
is  continuously  differentiable  in  x , there  is  a 

purely  imaginary  characteristic rooj.  A = iv(),v^  > 0, 

fo r a = a and  n o ot her  c h a_r a c t er ij^£i £ root 

A.  / A p , X o for  a = °q  satisfies^  A.  = m A ^ , m an 

integer ; 

(II.,)  r_  a ri)r  K > 0 ; : (-^,0]  -►  ]Rn  uniformly  continuous 

together  with  its  first  derivative  bounded  by  K, 
the  function  L(ct)d  + f(a,<j>)  lias  a derivative  in  u 
which  is  c o ntinuous  i n a , <J> . 

(1I3 ) 5(aQ)  t 0. 

Die  number  ?(a())  is  defined  in  the  following  manner. 
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Aq  - ' d>\)  " -iVf)  are  characteristic,  roots  of  liquation  (4.3) 
for  " n ,y  There  are  characteristic  roots  A(u)>Xfa)  of 
Equation  (4.3)  which  are  continuously  differentiable  in  ot , 

X (u()  l i v , and  two  Linear  independent  solutions  h exp  X(u)t, 
U exp  ^|x)t  of  Equation  (4.3).  Defining 


" C i ♦ ‘h  2 ^ = lic  >’  exp  X 

4.  ?(0)  --  I in  h exp  X(«)6, 

then  we  may  assume  (by  o change  of  basis  if  necessary)  that 


V0)  = *a(0)e 


11  (ot  d 9 


- «•  < 6 < 0 


lUat)  v03n  + B1(a) 


0 1 
1 0 


B,  (<*) 


Under  hypotheses  (IL)-(H^)  the  Bifurcation  Theorem  2.1  holds. 

The  proof  follows  the  one  given  in  [16]  under  less  general 
hypotheses.  The  appropriate  generalization  of  Theorem  2.2 
is  also  valid. 

In  some  applications,  it  may  be  more  appropriate  to  have 
other  spaces  of  initial  data.  For  any  p > 1,  one  can  consider 
Bp  defined  in  Relation  (4.2)  with  2 replaced  by  p (see  [16]). 
Another  interesting  space  is  the  following.  For  any  real 
number  y,  let 
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I 


k 


(4.9) 


— V Q 

0 ~ { 4s : ( -oc » 0 ] -►  Ik  continuous,  c p ( 6 ) a limit 

as  0 ->•  -oo} 


with  | T | y = sup,j{  |'H0 ) |oxp  - yO  } . If  V > 0,  the  essential, 
spectrum  of  the  semigroup  T(t)  has  the  same  propci ties  as  in 
B-,  (see  [10]).  The  Hop  I'  bifurcation  should,  therefore,  be 
valid  for  11  , y > 0.  A recent  preprint  of  Naito  ]70]  in- 
dicates it  will  be  possible  to  supply  a proof  of  this  result  in 
a manner  similar  to  the  one  given  in  Section  2.  Also,  this 
preprint  suggests  that  the  same  will  be  true  for  the  more 
general  spaces  of  initial  data  considered  hv  Hale  and  Kato  [10]. 

Considering  the  example  in  [16],  suppose  k is  nonnegative, 


con  t i mious  on  ( , 0 ] , 


and  consider  the  scalar 


J - f o 


equa  t i on 


(4.10) 


x(t) 


-a 


- r 

k(0  + r)f  (x(t+0))d0 

-ro 


where  r > 0,  a > 0 are  constants,  f is  continuous  together 
with  derivatives  up  through  order  two,  f(0)  = 0,  f(0)  = 1. 

We  wish  to  study  the  Hopf  bifurcation  for  liquation  (4.9)  as 
a fund  i on  of  (ot  ,r)  . 

In  [12],  Israelson  and  Johnson  and  in  [13],  Johnson  and 
Karlsson  proposed  the  liquation  (4.10)  as  a mathematical  model 

to  study  the  rhythmic  movement  of  some  plants.  experimental 
evidence  indicates  the  kernel  k should  have  the  form  shown 


I 


m 
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in  Figure  6.  'I  he  particular  form  of  k is  not  important 
be] ow . 


Figure  6 

Tho  first  problem  is  to  choose  the  space  of  initial 
functions  so  that  the  previous  theory  applies.  If  there  is 
a p > 0 and  constant  K such  that  |kf'3)|  < K exp  pO  , 
for  ■«  < 0 <_  0,  then  v/e  take  g(9)  = exp  uB , 0 < u < p and 


define 
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B = ft:  ( - o:* , 0 ] -»  IP:  t is  locally  integrable. 


0 


|t  | 4 = |4>  (0)  + I ey°|t(9)  | 2rl0  < <»} 

j - cti 


In  order  tor  the  right  hand  side  of  liquation  (4. 10) 


del 


(4.11)  F(r,o.,t)  USJ  -a 


- 1' 


1 .... 


k(9  + r)f  (t  («))d0 


to  he  continuous  in  (r,o,t)  and  continuously  differentiable 
in  u,t,  we  suppose  f satisfies 


(x)  | < a + hx  , I f ’ (x)  | < a1  + b1 


for  all  x and  some  constants  a,b,a^,bj  (see  [16]  or  [26]). 


The  function  F(r,a,t)  is  not  continuously  differentiable  in  r 
but  it  does  satisfy  Hypothesis  (H^)  above. 

The  linearization  of  liquation  (4.9)  around  x = 0 is 


(4.12)  x ( t ) = -a  k(0+r)x(t+e)de. 


-r 


The  characteristic  equation  is 


(4.1?)  X + a k(0+r)e  dO  = 0 


- r 


A6 


•Stech  [2?]  has  shown  that  the  liquation  (4.1?)  has  all 


j 


roots  with  negative  real  parts  if 


(4.14) 


-41 


Qk(e)cie  > 


Therefore,  no  bifurcation  can  occur  if  Relation  (4.14)  is 
satisfied. 

If  r > 0 , A = r A , o = r<x , then  liquation  (4.J2)  is 
equivalent  to 


(4.15) 


A + ac 


-A 


1 \ 9 

e k(r9)d(r0)  - 0. 


Stech  )23]  has  shown  that  for  « > v/2  there  is  an  r(J  > 0 
such  that  Equation  (4.15)  has  roots  with  positive  real  parts. 
This  certainly  suggests  there  must  be  a Uopf  bifurcation  for 
Equation  (4.19)  for  some  value  of  'i,r.  To  show  that  a smooth 
bifurcation  occurs  at  some  point,  one  must  be  able  to  say 
something  about  the  derivatives  of  the  roots  as  they  cross 
the  imaginary  axis,  generally  a very  difficult  problem. 

Lima  [16]  has  discussed  some  particular  coses  for  the 
function  k to  show  that  smooth  bifurcations  do  occur.  In 

the  class  of  all  continuous  kernels  k one  can  say  without 
any  computations  that  generically  there  will  always  be 
smooth  bifurcations. 

Of  course,  once  the  linear  equation  (4.13)  has  been 
analyzed,  the  existence  of  a Ilopf  bifurcation  will  also  be 
true  for  the  equation 
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as  well  as  more  general  equations. 

I r is  possible  also  to  consider  the  linear  equation  to 
have  the  form 


(4.16) 


x(t)  * -a  x(trO)dICft)  + r) 


where  K(d)  is  a nondecrens  ing  function  on  ( - ” , 0 ] . This 
allows  discrete  delays  as  well  as  integral  dependence  over 
the  past  history.  If  we  suppose  there  is  a p > 0 such 
tha  t 


e 1J  ^ d K.  f 0 1 < or' 


then  the  previous  theory  will  be  applicable  by  choosing  the 
space  11  . Of  course,  there  must  be  a few  technical  changes 
to  take  care  of  the  fact  that  <j>(-l)  for  example  is  not  a 
function  on  intograble  functions  on  (-«,()].  A more  convenient 
space  would  probably  Ire  C defined  in  Relation  (4.9). 

lire  results  of  Stech  [23]  are  also  valid  Tor  liquation  (4.16) 
Therefore,  bifurcation  problems  can  be  considered  for  either 


- r 

a ( t ) = -a  f(x(t  + 0))dK(6  + r) 

J ~ 00 


x ( t ) - -a  xft.  + 9)clK(9  + r)  (l  + x(t)] 


or  more  general  equations. 
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5.  Stabilizing:  effect  of  delays . . In  many  applications  of 
delay  equations,  an  increase  in  the  delay  tends  to  have  a 

destabilizing  effect  (see,  for  example.  Cur-lung  16),  May  [18]). 
In  population  dynamics,  it  lias  been  noted  that  it  can  have  a 
stabilizing  effect.  The  rate  of  growth  away  from  an  unstable 
equilibrium  can  lie  decreased  by  an  increase  in  the  delay  (see, 
for  example,  Cushing  [(>,  p.  35ff]).  lor  some  population  models 
representing  two  species,  Cushing  f 6 , p.  80ff]  proves  that  an 
unstable  equilibrium  cannot,  be  made  stable  by  increasing  the 
delay  although  the  instability  can  be  weakened  in  the  above 
sense.  He  alsc  gives  an  example  [6,  p.  99ff]  in  three 
dimensions  for  which  an  unstable  equilibrium  can  be  made  stable 
by  increasing  the  delays.  This  means  stable  Hopf  bifurcation 
can  occur  by  decreasing  the  delays. 

It  is  possible  also  to  give  an  example  in  two  dimensions 
where  unstable  equilibrium  can  be  made  stable  by  increasing  the 
delays.  The  purpose  of  this  section  is  to  give  an  example  of 

this  type. 

Consider  the  system 


(5.1) 


y(t)  = 

-CjjYU) 

- c12i(t)  * 

b^y  (t 

■v 

i(t)  - 

cny(tl 

- c22z(t)  ♦ 

b2y(t 

'r2} 

, > o, 

r.  > 0 

are  constants 

. No 

claim  is 

ij  J J 

made  to  the  physical  importance  of  this  equation,  but  it  is 
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similar  to  the  symbiotic  case  in  population  dynamics  (see 
Key  f i t z [14,  p . 2 7.V| ) . 


If  C = (c.j),  we  suppose 


(5.2) 


( c 1 1 + c22}  ' 4 det  C < °- 


This  implies  the  eigenvalues  of  C arc  - a ^ ± in  where 
a^  > 0 , a ? > (i;  that  is,  the  eigenvalues  of  C have  negative 
real  parts  and  are  complex  conjugate.  To  simplify  computations, 
also  suppose 


bl  B b2  = b’  rl  = r2  " 1 


By  a change  of  coordinates 


we  may  assume  System  (5.1)  has  the  form 


x(t)  = Ax(t)  + hx(t-r) 


(5.3) 


*al  " a2 


a 2 *al 


If  -a^  + b > 0 , r = 0 , the  ordinary  differential  (5.5) 
has  eigenvalues  with  real  parts  positive  and,  thus,  x = 0 is 


unstable.  It  ha?  been  shown  by  Tsen  [2?]  that  for  :i = ti / 2 
there  i «:  an  r > 0 such  that  the  zero  solution  of  liquation 

(5.3)  is  asymptotically  stable  for  r > r ^ . IV  e reproduce  some 
of  these  computations  hare. 

By  another  change  of  coordinates,  one  can  reduce  Equation 

(5.3)  to  a scalar  equation  in  a single  complex  variable  w, 

(5.4)  w = (-aj+ia.)w  '•  bw(t-r), 
for  which  the  characteristic  equation  is 

(5.5)  A = -a j + ia7  + be 


If  A = 0.  + iB}  a,  6,  real,  then 


(5.6) 


, -Cl  r 

01  = fll  + )e  C0S 

“ Ci  J* 

6 = a 2 - be  sin  Br 


or 


(5.7) 


0 _ „ . ru2  -2«r  2.  1/2 

B = a.,  ± [b  e - (0l  + a1)  ] 


a = -a^  + be  trcos  Be. 


Any  real  solution  of  these  equations  must  satisfy  G < a < ct  (r) 
where  Wg(r)  i s t^ie  uni(llJe  positive  zero  of 
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(5.8)  f(a,r)  = b2e'2Ur  - (o^)2. 

For  0 < « < ot^(r),  f(a,r)  > 0 and  f(u,r)r  1 0 ns  r >■  ® 
uniformly  in  a.  Choose  e > 0 such  that  | h | cos (V 2+ £)  < 
and  choose  r^  > 0 such  that  |f(a,r)r|  < e for 
0 < a < otp(r),  r > rQ  . Define 

B = a2  + [lr  e "r  - («+a1 ) “ ] i/ 

Since  a^r  = tt/2,  we  have  | Br  - v/2|  < e and 

a = -a j + be  cos  pr 
< -a j (1 -e  a) 

for  r > i'p.  This  latter  inequality  implies  ct  < 0 since 
a^  > 0.  This  proves  the  assertion  that  increasing  the  delay 
stabilizes  Equation  (5.1)  for  r^  - = r,  a^r  = ti/2,  b^  = b2- 

As  remarked  earlier,  this  implies  that  Equation  (5.1)  under 
the  above  hypotheses  and  when  subjected  to  nonlinear  perturba- 
tions which  vanish  together  with  their  first  derivatives,  will 
generally  have  a stable  Hopf  bifurcation  occur  as  the  delay 
decreases.  We  start  with  r > r^  where  the  zero  solution  is 
stable  and  decrease  r. 


6.  Nonautonomous  equations . The  previous  discuss  ion  is  con- 
cerned only  with  the  case  when  the  vector  field  in  Hie 
equations  is  independent  of  time  - the  autonomous  case.  When 
the  equation  is  nonautonomous,  the  problem  is  much  more 
difficult.  The  basic  difficulty  is  not  due  to  the  delays,  but 
is  a consequence  of  the  fact  that  it  is  difficult  to  determine 
the  effect  of  nonlinearities  in  resonance  phenomena.  The  best 
procedure  available  for  discussing  this  problem  locally  is 
general  transformation  theory,  averaging  methods  and  the  method 
of  integral  manifolds.  for  ordinary  differential  equations,  the 
basic  ideas  in  the  transformation  theory  go  back  as  far  as 
Liapunov  {151  (and  perhaps  further)  when  lie  was  concerned  with 
the  problem  of  determining  conditions  for  the  stability  of  an 
equilibrium  point  for  a nonlinear  system  when  the  coefficient 
matrix  of  the  linear  part  has  zero  roots  or  purely  imaginary 
roots.  Certain  aspects  of  the  method  of  averaging  were  also 
encountered  by  Liapunov.  The  more  general  theory  of  averaging 
as  well  as  the  theory  of  integral  manifolds  was  discovered 
by  Krylov  and  Bogoliubov  in  the  1930's  (see  Bogoliubov  and 
Mitropolski  [2]  or  Hale  [81).  The  transformation  theory  is 
extensive  (see,  for  example,  Malkin  (17),  Brjuno  [3]). 

The  generalization  of  these  ideas  to  delay  equations  of 
retarded  and  neutral  type  began  in  the  late  1950's.  An  idea 
of  the  contributions  in  this  development  can  be  found  in  [17] . 

At  this  time,  we  can  safely  say  that  most  of  the  results  that 
are  known  in  this  area  for  ordinary  differential  equations  are 
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known  for  delay  equations  or  the  necessary  machinery  is 
available  to  obtain  them. 

It  is  not  feasible  to  present,  the  elements  of  this  theory 
in  a few  pages  - a book  would  be  more  reasonable.  Therefore, 
we  are  not  going  to  discuss  delay  equations.  Fur thermore , we 
discuss  only  one  problem  for  a very  simple  ordinary  differ- 
ential equation.  To  avoid  any  technicalities,  we  also  assume 
all  functions  are  very  smooth  in  all  variables. 

Consider  the  second  order  equation 

(6.1)  x = Bx  + f(t,x,e)  . 

where  x = (x1  ,xj  elR",  f = (f^,f^,)  clR2,  B is  a 2x2 
constant  matrix,  the  real  parts  of  the  eigenvalues  of  B are 
zero , 


(6.2) 


a > 0 


the  function  £(t,x,e)  .is  T-periodic  in  t, 


(6.3) 


f ( t , 0 , e ) = 0 

f ( t , x , 0 ) = 0 ( | x | 2 ) as  I x I - n. 


Hypothesis  (6.3)  implies  x = 0 is  a solution  of 
Equation  (6.1).  Our  objective  is  to  discuss  the  behavior  of 


r SB 
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the  solutions  of  Equation  (6.1)  near  x ~ 0.  In  the  previous 
sections,  we  have  considered  this  problem  .for  f(t,x,e)  in- 
dependent of  t and  observed  that  it  was  possible  lor  a llopf 

bifurcation  to  occur  at  e = 0.  The  bifurcating  periodic 

2 

orbit  gave  a closed  curve  F in  1R  and  a cylinder  r x I’. 

n 

of  solutions  in  (x,t)-space  in  1R,‘  x 1R.  The  cylinder  is  obtained 
because  to  any  nonconstant,  periodic  solution  of  an  autonomous 
equation,  one  obtains  another  solution  by  a phase  shift 
(see  Fiaure  7) 
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If  the  bifurcating  periodic  orbit  is  asymptotically  stable, 

2 

then  the  cylinder  r * P.  is  asymptotically  stable  in  1R ' * IR. 

If  the  autonomous  differential  equation  is  subjected  to  small 

periodic  disturbances,  we  do  not  expect  to  preserve  the  same 

kind  of  orbit  structure  as  before  (.families  of  periodic 

solutions,  for  example),  but  we  might  expect  to  preserve  a 

surface  similar  to  the  cyl  inder  f >;  IK  which  would  be 

asymptotically  stable  and  completely  filled  up  with  solutions 

of  the  equation.  A surface  with  this  latter  property  is 

called  an  integral  manifold. 

Kith  the  above  motivation,  let  us  take  as  our  objective 

the  determination  of  integral  manifolds  in  a neighborhood  of 

the  set  {(x,t):  x = 0,  t £ IR)  which  are  like  cylinders.  More 

2 

specifically,  we  wish  to  find  a function  g:  IR  y-  IP  *•  JR  , g(t,ti) 
periodic  in  0 of  period  2n/o,  periodic  in  t.  of  period  T 
sucli  that  the  set 

S = { ( x , t ) c IR2  x IR:  x = g(t,e),  C t , 6 ) e IR2} 

is  an  integral  manifold;  that  is  if  (C,o)  c S,  then  the  solution 
x(t)  of  Equation  6.1)  through  (5 , o)  exists  for  all  t e IR 
and  (x(t),t)  e S for  t r IR. 

To  obtain  such  integral  manifolds,  let 


x 


1 


X 


2 


(6.4) 


= p sin  o9 , 


p cos  at' 
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to  oh  ta  i n 


(6  . S') 


» = 1 + = ( t , 0 , r , f ■ ) 
p = Rft,e,p,t) 


where 

5 ( t , 0 ,P,r  ) = (po)  ‘ ’ [ l^c.os  o0-  f 2^3.11  eO] 

(6.0) 

R(t,0,p,e)  = f,sj.n  00+  f,cos  00 

and  the  arguments  of  f .j  , f 2 lire  with  x given  in 

Relation  (6.4).  The  functions  5,R  satisfy 

5 ^ 0 ( ( 0 I +•  I e I ) 

(6.7)  R ( t , 0 , e ) - 0 

R(t , P,  0)  = 0(1  p|2), 

are  T-period  in  t,  2Tr/o-periodic  in  6. 

We  do  not  attempt  to  analyze  liquation  (6.5)  completely 
but  only  give  a special  result  to  illustrate  some  of  the 
basic  ideas.  Let 

r(t,6,P,e)  = 1 + n^t.O.e)  + ci2(tf0,e)P  + a?(t,0 ,e)p2  + 0( | p| ' ) 


(6.  it) 

R(t,0,P,c)  = 81(t,0,e)P  + e2(t,0,e)p2  + S3U, 0,<Qp3  + 0(|p|4). 
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Thc  objective  is  to  transform  to  now  coordinates 


(6.9) 


i - # + U0(t,0,E)  + U](t.,0,G)p  + U2ft,e,f.)p/‘  + o ( j P 1 ■5) 
r - p + v1(t,0,t)p  + v2(t,0,c.)p7  + 0(|p|3) 


where  all  functions  are  '(-periodic  in  t and  2 v/c  periodic 
in  t in  such  n way  that  the  new  equations  in  c. , i have  the 
f o i in 


(6.1 0) 


4=1+  VjCO  + U2(e)r  •*  u^(e)r2  + 0(|r|‘’’) 
r = v^fclr  + v2  ( e)  r 2 + v^(e)r5  + 0(|r|4) 


If  one  formally  substitutes  Relations  (6.9)  in  liquation  (6.5) 
and  equates  coefficients  of  r,  then  one  observes  that  each 
of  the  functions  uj,vj  satisfy  a partial  differential 
equation  of  the  type 


(6.11) 


If  * W - hCt.O) 


where  h(t,6)  is  T-perioclic  in  t and  2«/o  periodic  in  0 


To  obtain  a solution  of  liquation  (6.11)  which  is  T-periodic 

in  t and  2'l'/o  periodic  in  6,  and  must  impose  some  restrictions 

on  where  w = 2r/T,  In  fact,  if 
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u(t,0)  - l u^e 


Mt.o)  = l hkJle 


i (kmt'-AoO ) 


.i  (kwt  + P.aS ) 


then  we  can  uniquely  determine  u^  as 


uki  " '‘(k  *'■<>)  l'\f 


ht  - o 


is  k'i>  + 9.0  /■  0 for  every  h^  that  is  not  zero.  If 

we  take  u - 0.  Therefore,  if  we  assume  kw  + i.e  ^ 0 for 

kA 

those  k t 0,  9 f 0 which  occur  in  the  Fourier  series  for  h,  the 
function  u will  be  uniquely  determined  if 


2 rr  / o 2n/(o 


(6.12) 


J 0 


h(t,e)dtd0  = o 


In  general,  Relation  (6.12)  will  not  be  satisfied.  The  co- 
efficients y ..  ( e)  , v^.(e)  are  chosen  so  that  h will  satisfy 

Relation  (6.12).  Thus,  y.j(e),v. (e)  are  the  mean  values  of 
some  functions  explicitly  obtained  in  tlic  process  of  computing 
the  u.  and  v.  . For  more  details  on  methods  of  computation, 
the  reader  is  referred  to  [21],  [51],  [17]. 

If  we  assume  tire  nonresonance  conditions  are  satisfied 
for  m,ct  and  Equation  (6.5)  has  been  transformed  to 
Equation  (6.10),  then  we  can  state  the  following  result  which 
is  a consequence  of  the  theory  of  integral  manifolds. 
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Theo rein  (j_.  4 . I_f 

v'(0)  i 0 

(6.17)  v2(0)  - v’(0)  = 0 

v (0)  f-  0 , V 3 ( 0 ) v i (°)  < (l 

then  there  is  an  f()  > 0 ami  a function  yA1  , 0 , l1 ) , 0 < u < /t(), 

1/2 

T-per  iodic  _in  t,  2^/a- per  iodic  in  o> , g (t  ,6 , 0)  = [ >^(0)  (0)  ] 

such  that  the  set 

(6.14)  Se  = {(x,t)  e 1R2  x 1R;  x]  - /E  g ( t , 6 , /E) s in  06  , 

x9  = /e  g(t  ,0  , /e)cos  ct0  , 6 e 1H,  t.  e 1R } 

0 < e < is  an  integral  manifold  of_Fquation  (6.1).  Further - 

more,  this  integral  manifold  is  asymptotically  stable  if 
v- (0)  > 0 and  unstable  if  vj  (0)  < 0. 

If  v^(0)v|(0)  > 0,  the  same  conclusions  hold  on  existence 
if  -0^  < e _<  0 and  S£  is  unstable  if  vj  (0)  < 0 and 
asymp  tot. i cal  1 y_  stable  if  v j ( 0 ) > 0 . 

r, o t us  make  a few  more  remarks  about  second  order  ordinary 
differential  equations.  Consider  the  equation 

(6.15)  y = By  + F(t,y,e) 


A 


==— ’ — 1 
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l’(t,0,0)  = 0 
F(t,y,0)  = 0(|y|2). 


Assume  the  linear  equation 

)'  = By  4 h(t) 

for  a T periodic  function  i>  has  a unique  T- per  iodic  solution. 

Then  the  nonlinear  Equation  (6.16)  will  hare  a unique  T- 

£ 

periodic  solution  y (t,c)  in  a neighborhood  of  x = 0,  e = 0, 
y*  f t , 0 ) = 0 (see,  for  example,  [8]).  If 


y = y*(t,e)  + x 

then  x satisfies  an  equation  of  the  same  form  as  Equation 
(6.1).  Thus,  the  above  result  applies  to  Equation  (6.16). 

For  equations  with  delays,  results  similar  to  Theorem  6.1 
can  be  proved.  However,  much  more  machinery  is  needed.  One 
must  use  the  decomposition  theory  of  linear  autonomous  equations 
(see  Hale  [7,  Ch.  7,  9]),  averaging  methods  (see  Chow  and 
Mallct-Paret  [5]  or  Hausrath  [11])  and  the  generalization  of 
the  methods  of  integral  manifolds  to  delay  equations  (see,  Cor 
example,  Hale  [7,  Ch.  10,  Sec.  10.4]  for  some  references).  If 
one  is  very  familiar  with  these  methods  in  ordinary  differential 
equations,  the  ideas  used  in  delay  equations  are  not.  too 
surprising.  Before  tackling  these  problems  for  delay  equations, 


t 
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I would  suggest  tliat.  the  trader  he  able  to  supply  all  details 
in  the  previous  example. 
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